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The statistical model proposed in an accompanying paper is generalized to treat multiport scat- 
tering problems. Attention is first focused on two-port lossless systems and the model is shown to 
be consistent with Random Matrix Theory. The predictions are then tested by direct numerical 
simulation for a specific two-port cavity. Formulae are derived for the average transmission and 
reflection coefficients in terms of the port radiation impedance. The cases of cavities with multiple 
ports, and with a single port and distributed losses are compared. 
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I. INTRODUCTION 

In an accompanying paper a general statistical 
model is proposed to describe scattering of high fre- 
quency electromagnetic waves from irregular cavities 
(Here by high frequency we mean that the wavelength 
is substantially less than an appropriate characteristic 
length of the cavity.). This problem arises when one con- 
siders the coupling of electromagnetic energy into and 
out of complicated enclosures. The statistical approach 
is warranted when the exact details of the configuration 
are unknown or are too complicated to simulate accu- 
rately. The model proposed in Ref. fH gives a method 
for calculating the statistical distributions of important 
quantities that depend only on a small number of system- 
specific parameters. These system-specific parameters 
are the average spacing between resonant frequencies of 
the enclosure, the average quality factor, and the prop- 
erties of the coupling ports. It was shown Q that the 
key property needed to characterize the coupling port 
is its radiation impedance, that is, the impedance that 
would be seen at the port if the other boundaries of the 
enclosure were perfectly absorbing or were removed to 
infinity. This impedance can be used to "normalize" the 
fluctuating impedance at the port of the actual enclosure. 
With this normalization the fluctuating impedance has a 
universal (i.e., system-specific independent) distribution. 

The treatment in Ref. 0] focused on the case of a single 
port. The purpose of the present paper is to generalize 
these results to the case of multiple ports. A multiple 
port description is necessary if one is to describe the 
transmission of energy entering the enclosure to some 
other point in the enclosure. We will show that the mul- 
tiple port case can be treated using the same approach 
as the single port case. The main difference being that 
the multiple port description involves the introduction of 
matrices describing the coupling to the system. 
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Scattering can be described by either the impedance 
matrix, which relates voltages and currents at the ports, 
or by the scattering matrix, which relates the amplitudes 
of incoming and outgoing waves. In our model we fo- 
cus initially on the impedance matrix because its proper- 
ties can be related directly to the radiation impedances 
of the ports. Then based on our understanding of the 
impedance matrix, conclusions about the scattering ma- 
trix are drawn. 

The statistical treatment of wave scattering in com- 
plicated systems has been developed extensively in the 
physics community in connection with nuclear scattering 
1^] and scattering in mesoscopic systems such as quantum 
dots 0] and disordered conductors 3 ■ Many of the con- 
cepts developed in consideration of these problems have 
also been applied to electromagnetic scattering in mi- 
crowave cavities IS H B • The underlying approach 
is a description of the system based on random matrix 
theory (RMT) IQ. .11, .12]. Here, the specific systems un- 
der consideration are replaced or modelled by random 
matrices. These are matrices whose elements are inde- 
pendent Gaussian random numbers drawn from specific 
ensemble distributions. Two different ensembles that are 
considered are relevant to our problem. The Gaussian 
Orthogonal Ensemble (GOE) gives rise to a real symmet- 
ric matrix, where the variance of the diagonal elements 
is twice that of the off-diagonal elements. This applies 
to wave systems with time reversal symmetry (TRS), as 
would apply in media with real symmetric permittivity 
and permeability tensors. The Gaussian Unitary Ensem- 
ble (GUE) gives rise to complex, self-adjoint matrices, 
where the off-diagonal elements have independent real 
and imaginary parts. This applies to wave systems that 
have time reversal symmetry broken (TRSB). For elec- 
tromagnetic wave systems, this case is realized when a 
nonreciprocal medium such as magnetized ferrite or a 
magnetized cold plasma is present in the system. 

Random matrix theory, in spite of its high level of ab- 
straction, has been remarkably successful in predicting 
the universal statistical properties of wave systems. This 
includes the description of the distribution of resonances 
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in closed systems as well as the properties of scattering 
from open systems. One issue, that we address here is 
how to simultaneously account for the universal proper- 
ties as predicted by RMT and the system-specific prop- 
erties that depend on details of the coupling. In our 
approach this connection is made through the radiation 
impedance of the ports Id (T his connection has also been 
made by Warne et al. [13 )■ An alternative approach 
[T^ in the physics community, known as "the Poisson 
Kernel" , is based directly on the scattering matrix. This 
approach describes the distribution of fluctuations in the 
scattering matrix elements in terms of their average val- 
ues. Recently, Kuhl et. al.0 applied the Poisson Ker- 
nel approach to analyse the statistics of measured values 
of the reflection coefficeint for a microwave cavity. The 
measured results were found to be in agreement with pre- 
dictions based on the Poisson kernel. 

The organization of this paper is as follows. Section II 
generalizes the description of the quasi-two dimensional 
cavity introduced in Ref. [J] to describe the presence of 
multiple ports. Our treatment will be sufficiently gen- 
eral so as to treat the case of both isolated ports and 
non-isolated ports. We then discuss the changes that 
must be made when non-reciprocal elements such as a 
magnetized ferrite are present. In Sec. Ill we examine 
the statistical properties of the proposed impedance ma- 
trix, focusing on the two-port case. We also compare 
the model predictions with the results of computational 
electromagnetic calculations using HFSS. In Section IV 
we consider properties of the scattering matrix when the 
port coupling is imperfect. We find that the average value 
of the reflection coefficient is a function of the radiation 
reflection coefficient. Also, we consider the case of mul- 
tiple ports and verify the equivalence of distributed and 
diffractive losses in the case of large number of ports. 
Section V contains a summary of our results. Except for 
Section IV(B) all our considerations will be in the context 
of lossless cavities. 



II. GENERALIZATION OF THE MODEL 

A. The Impedance in the TRS Case 

Following Ref. 13] we consider a quasi-two dimensional 
cavity in which only the lowest order transverse magnetic 
modes are excited. The fields in the cavity are deter- 
mined by the spatially dependent phasor amplitude of 
the voltage VT{x,y). The voltage is excited by currents 
li drives at the various coupling ports, 

M 

(Vi + k')VT = -jfc/i?7o^M,/,. (1) 

1=1 

Here k = w/c, rjQ — -^/io/eo, h is the height of the cav- 
ity and an exponential time dependence exp(ja;<) has 
been assumed for all time dependent quantities. Each 
of the M ports is characterized by a profile function Ui 



centered at different locations and / dxdyui = 1. The 
phasor voltage at each port can be calculated as before, 
Vi = J dxdyuiVr =< UiVr > and is linearly related to 
the phasor currents Ij through the impedance matrix. 

To obtain an expression for the matrix Z, we expand 
Vt as before Q in the basis , the eigenfunctions of the 
closed cavity. The result is 

Z = -jkhvoj:P^^ (2) 

n ^ 

where the vector $n is [{uicj) ),...,{uM4'n)Y' ■ Us- 

ing the random eigenfunction hypothesis, we write (/>„ as 
a superposition of random plane waves. Thus the ele- 
ments of the Af-dimensional vector will be Gaussian 
random variables with zero mean. Elements of with 
different values of n corresponding to different eigenfunc- 
tions will be independent. However, for a given eigen- 
function the elements of may be correlated. This 
will be true, particularly, if two ports are close together, 
because the random superposition of plane waves leads 
to an autocorrelation function Jo(fc(5r) at two positions 
separated by 5r 16]. To treat correlations we write 

= L{k„)Wn, (3) 

where L is a non-random, as yet unspecified, M x M 
matrix that depends on the specific coupling geometry 
at the ports and may depend smoothly on fc„, and Wn 
is an M-dimensional Gaussian random vector with zero 
mean and covariance matrix equal to the M x M iden- 
tity matrix. That is we require that the components 
of the random vector Wn are statistically independent 
of each other, each with unit variance. Correlations be- 
tween ports are described by the off-diagonal components 
of L. The idea behind (0) is that the excitation of the 
ports by an eigenmode will depend on the port geometry 
and on the structure of the eigenmode in the vicinity of 
the ports. The dependence on the specific port geometry 
is not sensitive to small changes in the frequency or cav- 
ity configuration and is embodied in the matrix quantity 
L{k). The structure of the eigenmode in the vicinity of 
the ports, however, is very sensitive to the frequency and 
cavity configuration, and this motivates the statistical 
treatment via the random plane wave hypothesis. From 
the random plane wave hypothesis, the excitation of the 
port from a given plane wave component is random, and, 
since many such waves are superposed, by the central 
limit theorem, the port excitation is a Gaussian random 
variable, as reflected by the vector Wn- In Paper I, we 
have derived a result equivalent to Q for the case of a 
one-port with a specific model of the excitation at the 
port (namely, a vertical source current density Iu{x,y)z 
between the plates). Our derivation here will be more 
general in that it does not depend on a specific excita- 
tion or on the two-dimensional cavity configuration used 
in Paper I. Thus this derivation applies, for example, to 
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three dimensional cavities, and arbitrary port geometries. 
From ^ and ||2J) we have for the Z matrix 



jkhrjo ^ 



L{kn)Wn'W^L^ (fc„) 
1,2 jU2 



(4) 



We now take the continuum hmit of Ql and average 
over Wn, 



POO 

{Z) = ~j / khi^^L{k') 
Jo 



fc2 - (A:')' 



:L^{k')^. (5) 



where A is the averaged spacing in fc^ values. We note 
that the continuum limit is approached as the size of 
the cavity is made larger and larger, thus making the 
resonance spacing (fc^^^ ~ ^n) approach zero. Thus, the 
continuum limit corresponds to moving the lateral walls 
of the cavity to infinity. Using our previous one-port 
argument as a guide, we anticipate that, if the pole in 
Eq. lO at fc'2 = fc2 is interpreted in the causal sense 
(corresponding to outgoing waves in the case with the 
walls removed to infinity), then {Z) in (0 is the radiation 
impedance matrix. 



{Z) = Zn{k) = RR{k) + jXR{k), 



(6) 



where V — Z^ik)! with V the M-dimensional vector of 
port voltages corresponding to the M-dimensional vector 
of port currents /, in the case where the lateral walls have 
been removed to infinity. With the above interpretation 
of the pole, the real part of Eq. Q yields 



Rl^{k) = TTkhT]oL{k)L'^{k)/A. 

Thus, Eq. (0) becomes 



h.2 _ h.2 ■ 



(7) 



(8) 



where {wnW^) = 1m- (Note that the formula for A is 
different in two and three dimensions.) In the case of 
transmission line inputs that are far apart, e.g., of the 
order of the cavity size, then the off-diagonal elements 
of Zfi are small and can be neglected. On the other 
hand, this will not be the case if some of the transmis- 
sion line inputs are separated from each other by a short 
distance of the order of a wavelength. Similarly, if there 
is a waveguide input to the cavity where the waveguide 
has multiple propagating modes, then there will be com- 
ponents of V and / for each of these modes, and the 
corresponding off-diagonal elements of Zr for coupling 
between these modes will not be small. 

For the remainder of the paper, we will assume iden- 
tical transmission line inputs that are far enough apart 
that we may neglect the off-diagonal elements of Z^. As 
before, we will take the eigenvalues k"^ to have a distri- 
bution generated by RMT. Because the elements of Z 
depend on the eigenvalues fc^, there will be correlations 
among the elements. In the lossless case the elements of 



the Z matrix are imaginary, Z = jX, where X is a real 
symmetric matrix. Consequently X has real eigenvalues. 
We will show in Sec. Ill that the distribution for individ- 
ual eigenvalues of X is Lorentzian with mean and width 
determined by the corresponding radiation impedance. 



B. Effects of Time-Reversal Symmetry Breaking 
(TRSB) 

In the time-reversal symmetric system, the eigenfunc- 
tions of the cavity are real and correspond to superposi- 
tions of plane waves with equal amplitude waves propa- 
gating in opposite directions as in Eq. (13) of paper [Jl, 
which is recalled as follows 

— N 



^ N^^oo V A/V^^^E exp(jfc„ei • x + jOi)}, (9) 

i—l 

where a^, Oi and e*; are random variables. If a non- 
reciprocal element (such as a magnetized ferrite) is added 
to the cavity, then time reversal symmetry is broken 
(TRSB). As a consequence, the eigenfunctions become 
complex. Eq. @ is modified by removal of the operation 
of taking the real part, and the in Eq. (12) of pa- 

per m also become complex. In practice, there exists a 
crossover regime for the transition from situations where 
time reversal symmetry applies to those it is fully-broken. 
An interested reader might refer to situations where dis- 
cussion in Ref. |0| and the references therein. However, 
in this paper, we will discuss only the case when time- 
reversal symmetry is fully broken. In this case we find 



= [Ai?fl,(fcn)] 



1/2 



Win 



(10) 



where = (w^^^ -|-jw^^^)/\/2 and w^^-* and w^*] are real, 
independent Gaussian random variables with zero mean 
and unit variance. The extra factor of \/2 accounts for 
the change in the normalization factor in Eq. 0, re- 
quired when the eigenfunctions become complex. Fur- 
ther, transpose w^, in Eq. ||SJ| is now replaced by the 
conjugate transpose . 

A further consequence of TRSB is that the distribution 
of eigenvalues is changed. The main difference is the be- 
havior of P{s) for small s. In particular, the probability 
of small spacings in a TRSB system (P(s) ~ s^) is less 
than than of a TRS system (P(s) ~ s). 

For the sake of simplicity, we will assume all the trans- 
mission lines feeding the cavity ports are identical, and 
have the same radiation impedance, Z^ — Rji + jX^ — 
{Rr+ i X b)1 M , where Rn and Xn are real scalars. Anal- 
ogous to the one port case, we can define a model normal- 
ized reactance matrix = Xij/Rn for the case RR{kn) 
constant for n < and RR{kn) = for n > A^, 



1.^ 



fc2 



(11) 
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where P = P/A, win = {win +jWin)/V2, w'fj and 
are real, independent Gaussian random variables with 
zero mean and unit variance, E{w*j^Wjn) — Sij. Note 
that a unitary transformation, ^' = USJJ\ returns 1)11(1 
with Win and Wjn replaced by and where w'^ = 
Uwn- Since a unitary transformation does not change the 
covariance matrix, E{winW*n) = E[w[^w'*n) = Sij, the 
statistics of ^ and of ^' are the same; i.e., their statistical 
properties are invariant to unitary transformations. 



III. PROPERTIES OF THE IMPEDANCE 
MATRIX AND EIGENVALUE CORRELATIONS 
FOR LOSSLESS CAVITIES 

The universal fluctuation properties of the Z matrix 
can be described by the model matrix ^ij specified in 
Eq. (|ll|l . In the TRS case the Wjn are real Gaussian 
random variables with zero mean and unit width and the 
spacings satisfy Eq. (16) in Ref. Q. In the TRSB case 
the Wjn are complex and the spacings between adjacent 
kl satisfy Eq. (17) in Ref. Ilj. 

In the case under consideration of multiple identical 
ports, ^ij will have a diagonal mean part ^6ij for which 
all the diagonal values are equal. The eigenfunctions of 
^ij — £,5ij and of its fluctuating part will thus be 
the same. Consequently, we focus on the eigenvalues of 
the fluctuating part. 

We initially restrict our considerations to the two- 
port case. We recall that for the lossless one-port case 
there is no difference in the statistics of the normalized 
impedance ^ for the TRS and TRSB cases. In both cases 
, it is Lorentzian with unit width. In the lossless two-port 
case, however, essential differences are observed when 
time reversal is broken. Using Hll|) . we generate 10^ real- 
izations of the 2 by 2 matrix ^ in both the TRS and TRSB 
cases, again for N = 2000 and fc^ — 1000. In this test we 
generated spectra based on an independent spacings Q. 
For each realization we compute the eigenvalues of the ^ 
matrix. We find that in both the TRS and TRSB cases 
the eigenvalues of the matrix are Lorentzian distributed 
with unit width. That is, histograms of the eigenvalues 
generated according to the TRS and TRSB prescriptions 
are identical. However, if we consider the joint proba- 
bility density function (PDF) of the two eigenvalues for 
each realization, then differences between the TRS and 
TRSB cases emerge. We map the two eigenvalues ^i, 
j = 1 or 2, into the range [7r/2,7r/2] via the transfor- 
mation 9i = arctan(^i). Scatter plots of 62 and 9i for 
10^ random numerical realizations of the ^ matrix are 
shown in Fig. ^a) for the TRS case and in Fig. Cfb) for 
the TRSB case. The white diagonal band in both cases 
shows that the eigenvalues avoid each other (i.e., they are 
anti-correlated). This avoidance is greater in the TRSB 
case than in the TRS case. The correlation, 

corrie.M^^^^^4=Mm. (12) 




FIG. 1: (a) Scatter plot of Oi vs 62, in the TRS case, (b) 
Scatter plot of 6\ vs 62 in the TRSB case.(c) Scatter plot of 
6\ vs 62 from the HFSS simulation, with 100 realizations and 
sweeping frequency from 6.75GHz to 8.75GHz. 



is numerically determined to be -0.216 for the TRS case 
and -0.304 for the TRSB case. 

^From the construction of the ^ matrices for the TRS 
and TRSB cases their statistical properties are invariant 
under orthogonal and unitary transformations, respec- 
tively. Random matrix theory has been used to study 
these rotation-invariant ensembles and predicts the joint 
density function of Oi and O2 11] to be, 

Pp{e^,62)o,\e'^'^ -e'^'^f, (13) 

where /3 = 1 for the TRS case and /3 = 2 for the TRSB 
case. Note that based on Eq. H13|) . the probability density 
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function for one of the angles P{8i) = J rf^2-P(^i,^2) is 
uniform. From the definition 9 — arctan^, this is equiva- 
lent to the eigenvalues of the ^ matrix having Lorentzian 
distributions (P^(6) = PeiO^)\deJd^,\ = \dejdi,\/2n). 

The correlation coefficients calculated from the numer- 
ical results in Figs.^a) andQIb) are consistent with the 
predictions of the random matrix theory from Eq. (|13|l . 
that is, -0.216 for the TRS case and -0.304 for the TRSB 
case. This implies that the distribution of spacings and 
the long range correlations in the eigenvalues of the ran- 
dom matrix, which are ignored in the construction of the 
fc^ in the above computation are not important in de- 
scribing the statistics of lossless impedance matrices. As 
we have discussed in 0], these correlations could be in- 
cluded using a sequence of fc^^ generated by the eigenval- 
ues of a random matrix. (We note that 1] , lossy cavities 
yield statistics that are different in the TRS and TRSB 
cases.) 

Now we test these predictions for numerical simula- 
tions of the chaotic cavity considered in paper 0]. We 
use the HFSS software to calculate the cavity impedance 
matrix and radiation impedance matrix for a 2-port case. 
We locate the two ports, at the positions (x, y) = (14cm, 
7cm) and {x, 2/)=(27cm, 13.5cm). As in we also in- 
clude the 0.6 cm cylindrical perturbation which is located 
alternately at 100 random points in the cavity, and we 
numerically calculate the impedance matrix for 4000 fre- 
quencies in the range 6.75GHz to 8.75GHz. We obtain a 
normalized Z matrix, which is analogous to the ^ matrix 
defined in Eq. (|ll|l according to 

^hfss = R^\lm[Z,av] - I2XR), (14) 

where I2 is the 2 by 2 identity matrix, Zcav is the 2 by 2 
impedance matrix calculated by HFSS, and Xn and Rr 
are the radiation reactance and resistance for a single 
port. For each realization of ^hfss we calculate its eigen- 
values = tanOi, i = 1, 2, and plot the values on the Oi 
vs. 02 plane, as shown in Fig. ^c). The anti-correlation 
of the angles is seen in the figure, and corr{6i,92) from 
(|12ll is -0.205, which is comparable with what we expect 
for the TRS case, -0.216. 

So far we have focused on the eigenvalues of the 
impedance matrix. The eigenvectors of Z are best de- 
scribed in terms of the orthogonal matrix whose columns 
are the orthonormal eigenfunctions of Z. Specially, in 
the TRS case, since ^ is real and symmetric, 

where is the transpose of O, and O is an orthogonal 
matrix, which we express in the form 

O^fcosrj sinr;\ 

K—smr] cos rj J ^ ' 

A scatter plot representing the joint pdf of the angle r\ 
and one of the eigenvalue angles d\ is shown in Fig.|2Ial). 




FIG. 2: Scatter plot of 7? vs for (al) the model impedance in 
the TRS case, (bl) the model impedance in the TRSB case, 
and (cl) from the HFSS simulation. Plots (a2) and (a3) [(b2) 
and (b3), (c2) and (c3)] show conditional probability for Q and 
for 77 for the model TRS case [model TRSB case, the HFSS 
simulation] . 

In analogy to how we obtain the realizations used in 
Fig. 2 in [l|, this plot is obtained by inserting random 
choices for the k\ and win in 1)11(1 . Notice that we have 
restricted 77 in Fig.|2Ial) to the range < 77 < 7r/2. This 
can be justified as follows. The columns of the matrix 
O in (|16|l are the eigenvectors of We can always de- 
fine an eigenvector such that the diagonal components of 
O are real and positive. Further, since the eigenvectors 
are orthogonal, one of them will have a negative ratio 
for its two components. We pick this one to be the first 
column and hence this defines which of the two eigen- 
values is B\ . The scatter plots in Fig. show that the 
restriction on rj maintains the symmetry of 61 and 62, 
vis. Pp{6i,e2) = Pp{e2,0i). Also in the ¥ig.^a2) (and 
(a3)), we plot the conditional distribution of 9 (and if) for 
different values of ry (and 6*). As can be seen, these plots 
are consistent with 77 and 9 being independent. This is 
also a feature of the random matrix model ^2 • This in- 
dependence will be exploited later when the S matrix is 
considered. 

For TRSB systems, the ^ matrix is Hermitian = £,* ■ 
A unitary matrix of eigenvectors that diagonalizes it can 
be parameterized as 

U^i COST? sin 77e<\ 
sm -qe ^ cos rj J 

Thus, there is an extra parameter C, characterizing the 
complex eigenvectors of the ^ matrix. According to ran- 
dom matrix theory, the eigenfunctions and eigenvalues 
are independently distributed, i.e. 77 in the U matrix 
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should be independent of 0i, 02- This expectation is con- 
firmed in Fig. I3b) where a scatter plot of 9i vs rj and 
conditional distributions of 9 and of rj are shown. 

Again, we test the independence of 6 and 77 with 
HFSS calculations. Using the ^hfss matrix obtained from 
Eq. H14|l . the angles 6 and 77 can be recovered from the 
eigenvalues and the eigenvectors of the £,hfss- With the 
ensemble generated by sweeping the frequency from 6.75- 
8.75GHz and considering 100 different locations of our 
cylindrical perturber, we obtain the joint distribution of 9 
and 77 in Fig.EJcl) as well as their individual distributions 
in Fig.[2tc2) and (c3). Here we see that the distributions 
are qualitatively similar to those of the model impedance 
matrix in the TRS case. However, there are significant 
departures which need to be investigated. It is likely that 
these are the result of the same strong multi-path inter- 
ference which gave rise to the reactance variations in the 
one port case shown in Paper 0. 



IV. AVERAGED REFLECTION COEFFICIENT 
IN MULTI-PORT CASE 



In this section, we use our knowledge of the statis- 
tical properties of the Z matrix to deduce properties 
of the S matrix, particularly for the ensemble average 
of the reflection coefficient (1511 p). For a system with 
two ports, in the lossless case considered here we note 
{\Su\') = l-{\Sn\') 

According to the previous section, for the case of non- 
perfect coupling, model of the cavity impedance matrix 

"1/2 "1/2 " 

can be expressed as Z = i?^ £,Rr + jXr, where Zji is 
the 2x2 radiation impedance and ^ is a 2 x 2 random 
matrix generated according to Eq. pi|l . If the incom- 
ing frequency is restricted in a narrow range, the radia- 
tion impedance Zn is essentially constant. In this paper 
we assume that identical ports are connected to identi- 
cal transmission line, i.e., Zj^ and the transmission line 
characteristic impedance Zq are diagonal matrices with 
equal diagonal elements. Thus, we obtain the expression 
for the S matrix, S = {Z + Zo)~^(^ - Zq), 

s = [{iR^ + nxi2) + i2r'[hR^ + nxi2)-i2l (18) 

where = Rr/Zq, = ^R/Zi) are scalars and I2 is 
the 2x2 identity matrix . These two parameters, as we 
show later, fully specify the coupling effects on the wave 
transport process. The special case of perfect coupling 
corresponds to ^r = 1 and 7x = 0. 



A. Lossless Two-port Case 

We recall that for TRS systems the reactance matrix 
X is real and symmetric, and can be diagonalized by 
an orthogonal matrix O, Eq. Hlt)|) . If identical ports are 
connected to identical transmission lines of characteris- 
tic impedance Zq, then the scattering matrix S is also 



diagonalized by O, and we can write 



S = 



,301 









(19) 



The scattering phases <j)i and (j>2 are then related to the 
eigenvalue angles 9i by formulas analogous to the one- 
port case, tan(7r/2 — (j)i/2) — jr tan 9i + ■yx- 

Substituting Eq. H16|l for O in (|19|l and multiplying the 
matrices, we obtain 

l^iil^ = cos^?7 + sin''r7 + 2cos^77sin^77cos(0i— 02)- (20) 

We can now compute the expected value of the square 
of 1 5*11 1 by assuming that 77 is independent of the an- 
gles (t)i and 02 and is uniformly distributed, which yields 
(cos* 7] + sin* 7;) = 3/4, 2 (cos^ rj sin^ 77) = 1/4 and 



\Sn\') 



3+ l(eos(0i-02)). (21) 



Assuming the angles 9i and ^2 are distributed according 
to Eq. (|13|l and using the relation between (j)i 2 and ^1.2, 
evaluation of (cos(0i — ^2)) is carried out in Appendix. 
The result is 



(|5iiP> 



l-'-^-^^^M^ln'-M, (22) 



8\pr\ 



where the "the free space reflection coefficient" pR is de- 
fined as the same way in the Paper 0], 



PR 



IR + iix ~ 1 

iR + hx + 1 ' 



(23) 



We first check the asymptotic behavior for the power 
transmission coefficient T = 1 — \S\\\^ implied by the 
formula (|^ . In the non-coupled case, \pr\ = 1, i.e., 
all the incoming power is reflected, and we obtain from 
(|22|l (r) = 0. On the other hand, in the perfect coupling 
case, \pr\ = 0, ln[(l + \pr\)/{1 - \pr\)] in the ^ can 
be expanded as 2{\pr\ - \pRf' /?>). Therefore, (T) = 1/3. 
This is consistent with the result in Ref. (-R) = 2(T). 
That is, in the perfect coupling case the average of the 
reflected power is twice that of the transmitted. 

Eq. H22I) shows that the averaged power reflection and 
transmission coefficients only depend on the magnitude 
of PR, and not its phase. A plot of (15*11 p) versus \pr\ 
is shown in Fig. E^a). Also shown are data points ob- 
tained by taking 10^ realizations of the impedance matrix 
(|ll|l with eigenvalue statistics generated from TRS spec- 
trum and computing the average of |5iip for different 
combinations of ^r and 7x characterizing the radiation 
impedance. The data confirm that the average of |5iip 
depends only on the magnitude of the free space reflec- 
tion coefficient and not its phase. 

In the TRSB case, the eigenvalues of the X matrix are 
still real, but the eigenvectors are complex. In this case, 
Eq. H19|) is replaced by 



S = U 



,3<Pl 



C,0<P2 



(24) 
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where the unitary matrix U is given by Eq. (|17|l . Multi- 
plying the matrices in Eq. H24|l . we find the same expres- 
sion for IS'iip, Eq. if^ . as in the TRS case. The average 
of IS'iip will be different in the TRSB case because of 
the different statistics for 77, Oi and O2 which character- 
ize the eigenfunctions and eigenvalues of the impedance 
matrix. In particular, 77 has a distribution, arising from 
the SU(2) group 18], 



F„(77) = |sin(277)|, 



(25) 



which yields (cos** 77 -I- sin^ 77) = 2/3, 2(cos^ 77 sin^ 7;) — 
1/3, thus. 



\Sii?) = 



cos 



-02)). 



(26) 



Recalling that 9i and 62 are distributed according to l(T^ 
with /3 = 2, this results in a different set of integrals (see 
Appendix). The result is 



(1^1 



= 1 - 



6 



(27) 



which depends only on the magnitude of the free space re- 
flection coefficient. A plot of (ISnp) from Eq. (|27l) versus 
\pb\ is also shown in Fig.|3{a), along with data point ob- 
tained by taking 10^ realizations of the TRSB impedance 
matrix Hll|l generating from random numbers and com- 
puting the average of IS'np for different combinations 
of 7fl and characterizing the free space impedance. 
Once again, the data collapse to the curve predicted in 
Eq. 123 • 

We now test the relation between (jS'iip) and \pb\ 
with the impedance matrices we obtained from the HFSS 
two-port calculations. We can vary the transmission line 
impedance Zq and generate (jSup) and |p_r|. However, 
the range of \pb\ values accessible doing this is limited 
because of the large inductive radiation reactance associ- 
ated with the coupling port. To extend the range of \pr\ 
we add a shunt susceptancc Y — (jtoC) in parallel with 
each port. This results in a modified cavity impedance 
matrix Z^^^ = {Z^J^^ + jujCl2)~^ ■ We then form the 
scattering matrix 



Zq) 



(28) 



The corresponding free space reflection coefficient is gen- 
erated by Z'j^ = (Z^^ + jwC)~i and \pr\ = \Zj^ + 
Zq\~^\Zj^ — Zq\. By choosing appropriate combinations 
of ujC and Zq, we can achieve a range of \pb\ values be- 
tween and 1. For each \pfj\ value, we average \Sii\'^ 
over frequencies and realizations and plot the points on 
Fig. |2Jb) . These compare favorably with the theoretical 
result (solid curve) based on the random matrix theory 
results. 



B. M-port Case, M > 2 

Using the random coupling model and assuming 
perfect coupfing -fj^ = I, -fx = ^ (i.e. \pr\ = 0), we 







TRS data 
X TRSB data 
--- Eq. (24) 
Eq. (29) 









0.2 0.4 0.6 0.8 




0.85 
— 0.8 
0.76 
0.7 
0.65, 



FIG. 3: (a) Numerical simulation for the average reflection 
coefficient (ISnp) vs magnitude of pR defined in Eq. 1231 for 
the TRS and the TRSB system, taking lO'' realization of the 
impedance matrix, 30 uniformly spaced values of from 0.1 
to 3, and 31 equally spaced values of yx from to 3. (b) 
Average reflection coefficient (liSiip) vs \pr\ using the cavity 
impedance and radiation impedance from HFSS calculation 
and varying the values of Zo and the capacitive reactance Y 



have simulated the S matrix for cases of two to seven, 13 
and 57 ports. The results for the average reflection and 
transmission coefficients were found to satisfy: 



TRS 



M+l 



and 



TRSB 



« = J, 



« =J, 



(29) 



(30) 



where M is the number of ports connecting the cavity 
to transmission lines. It seems that, in the TRS case, 
the input waves "remember" their entry port and have 
a preference for reflection through it (this is related to 
the concept of "weak localization" reviewed in 0|). In 
contrast, for the TRSB case, the waves behave as if they 
forget through which port they entered the cavity, and 
thus all the ports have equal probability of being the 
output for the waves. 

It was shown by Brouwer and Beenakker [23| that scat- 
tering in multiport lossless systems can be related to that 
in a single-port, lossy system. It was proposed that the 
introduction of N' (N' >> 1) fictitious ports into the 
scattering matrix of a lossless system would give equiva- 
lent statistics for the refiection coefficient seen at a single 
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port as would be obtained for a single port model with 
a uniform internal loss. Considering a system with M 
ports all perfectly matched, we can pick port 1 as the in- 
put and consider the other ports as a form of dissipation. 
Due to the energy escaping from the other (M — 1) ports, 
we will obtain a reflection coefficient Sn with magnitude 
less than one, which is similar to that obtained in the 
one-port lossy case (i.e., with losses due to finite wall 
conductivity). The cavity impedance seen from port 1, 
Zi, is calculated from S'n, one of the elements from the 
M by M scattering matrix. 



Zi = Rr 



l + Sii 



jXr. 



(31) 



When normalized by the radiation impedance this cor- 
responds to a complex impedance Cm = (l + 'S'ii)/(l — 
•^ii) = P + jS.- On the other hand, we can generate the 
lossy one-port impedance ( from Eq. (48) in Ref. 0, 
modelling the lossy effect by adding a small imaginary 
term to the frequency 0. We can then compare the 
statistics of C from the lossy one port and Cm from the 
M-port lossless case (We note that approximate analytic 
formula for the distributions of the real and imaginary 
parts of C have recently been given by Fyodorov and 
Savin ^21]). An appropriate value of the damping pa- 
rameter in the one port case, k^a (cr = l/Q), can be 
determined so that the average value of jS'iip in the 
lossy case is equal to 2/ (M -I- 1) for the TRS case (or 
1 /M for the TRSB case) . Then we can compare the real 
and imaginary parts of the impedances obtained in the 
two different ways. In Fig. 4, we include the results for 
the three different number of ports, M—A, 13 and 57, 
and the corresponding one port result. For M = 4 we 
note that the distributions are similar but clearly not 
the same. However, for M—13 or 57, the distributions 
for C and Cm are much closer. Thus, we confirm that 
distributed damping and a large number of output chan- 
nels are equivalent so as to affect the distribution of the 
sub-unitary scattering matrix. 

We now briefly discuss the multiport case with M > 2 
and with mismatch (|p_r| > 0). As long as the assump- 
tion that the eigenfunctions (77) and the eigenvalues {9 
or (j)) are independent is still true, (IS'iip) is related to 
the mismatch only through (cos((/)fe — 0/)), similar to the 
expression in Eq. H21|l . The same series of steps spec- 
ified in the Appendix can be carried out to show that 
(cos(0fe - 4>i)), as well as depend only on |p_r| 

(and are independent of the phase of pr). We have ver- 
ified this by numerical simulation using the impedance 
matrix generated from Ullfl with up to seven channels. 



V. SUMMARY 

We have generalized our random coupling model pro- 
posed in Ref. 1] to the multiport case. A similar 
impedance normalization is applied to obtain the statis- 
tical properties of the multiport chaotic scattering ma- 
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- - 57 Ports 
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FIG. 4: Comparison between the impedance obtained from 
the one-port lossy case and the multiple lossless case, (a) for 
the real part of the impedance; (b)for the imaginary part of 
the impedance. 



trix. The correlation coefficients between eigenvalues are 
calculated explicitly and agrees with the random matrix 
theory. 

We further incorporate the coupling parameters 7^ 
and 7x into the formulation of multiport scattering ma- 
trices and present the formula for the averaged reflection 
coefficients versus different values of coupling strength. 
We find that \pr\, which is a function of the two pa- 
rameters above, characterizes the transport process. For 
different pairs of (7^, ^x), as long as they yield the same 
value of \pr\, the corresponding averaged reflection coef- 
ficients are the same. This observation may offer a useful 
criteria for cavity design. 

Using HFSS, we test the conclusions above using 
impedance data calculated from direct numerical solution 
of Maxwell Equations. The agreement between the nu- 
merical results and the theoretical predictions convinces 
us that our approach of impedance normalization suc- 
cessfully recovers the statistical ensemble for chaotic scat- 
tering in the multiple port case. 
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APPENDIX: EVALUATION OF {\Siif) 

In this appendix, we will start from the one-port case, 
and obtain an expression for the phase of S in term of the 
reflection coefficient p/j defined in Eq. (|23|l . Then, using 
Eq. (fT^ . we can evaluate (cos((/)i — ^2)) for the two-port 
in the TRS and TRSB cases. 

In the one-port case, S can be expressed as 

Z-Zo 



Z + Zo 

j{ix +hR) - 1 
j{ix +hR) + 1' 



(A.I) 



where ^ is a zero mean, unit width, Lorentzian random 
variable, which can be written as, 



^ = tan( 



(A.2) 



with 9 uniformly distributed in [— 7r/2, 7r/2]. Putting 
Eq. I|A.2(I into Eq. (|A.1|I . we get 

^ {iR + ilx - l)e^'^ - (7fl - ilx + l)e~^'^ 
{iR + ilx + l)e^" - (7i?, - j7x - l)e-J" ' 
Introducing such that 

IR + ilx - 1 = Pr{ir + j7x + 1), (A. 4) 
and defining 

IR + ilx + 1 

we obtain a compact expression for (j) in term of 9 and 



PR 



1 - pIjC 



-j(2e+a) 
/ 1 



\PR\e- 



(A.6) 



where 29' = {29+a+Tr+(j)B.)- Since a and 0/{ depend only 
on the coupling coefficient 7/j and , and 20 is uniformly 
distributed in [0, 27r], the angle 29' is also uniform in 
[0,27r]. Thus, 



P^{4>) = P2e'{29 
1 



,.,d{29'). 



(A.7) 



27rl + |pfi|2-2|pfl|cos((/> -</./!)■ 

The relation between (f) and 29' also holds true for 
multi-port cases. Furthermore, from the joint probability 
density function of 26*1 and 2^2 in Eq. , which is only 
a function of the difference of two angles, we find that 
29'i and 2^2 have the same joint distribution specified in 
Eq. H13|) . Thus we can evaluate 

(cos(<?!)i - <?!)2)) = Re[e^'l''-^'l'^] 



Re\ 



+ \pr\ 



+ \pr\ 



(A.8) 



by using the joint distribution of 29'i and 2^21 
Pp{29i,292) oc |eJ20; _ e^^e^p^ ^^lem (3 = 1 corresponds 
to the TRS case, and /? = 2 for TRSB case. 

Introducing i/ji = 29'^, '02 — 26*2, and their difference 
■0_ = "01 — '02 , we obtain for the average of cos(0i — 02), 



(COS(01 - 02)) = 



i?e[ 



i?e| 



2~-f^W'^2) 



{2nY 



\pr\ e^'^^ + lpRl 



l + \pR\e-i^^ l + \pR\e^^^' 



2tt 

2ix 



■P(^_) 



(A.9) 



i>2 e'-''('^-+^^) + \pr\ 
'0 27r 1-1- |pfl|e-J('/'-+'/'2) 



e^'^^ + \pr\ . 
1 + |p/?,|eJ''/'2^' 

The inner integral can be calculated by introducing a 
complex variable z = e-''^'^ in terms of which the inner 
integral becomes 



1 

27j 



dzf{z) 



unitcircle 



z{z+\pR\e~^^-)' 



(A.IO) 



where /(z) = (|pfl|z+e-J'^- )(z+|pflj)/(l-f z|pfl,|). Evalu- 
ating this integral via the residues at the two poles within 
the unit circle, z — and z = —\pR\e~^'^~ , we obtain 



{coa[(Pi - 02)) = / ^— 



[1- 



° . (A.ll) 

(l-|pKn(l~cosV>_) 

1 + |pfl|4-2|pfl|2cOsV^_^' 



For the T RS case, Pv-lV^-) = 7r| sin(-0_/2)|/2, and 
Eq. HA.lljl yields 



(COS(01 - 02)) = 



|p^|4 + 2|p^|2-l 



(A.12) 



Mpr[ 



In- 



1 - \pr\ 



For the TRS B cas e, Pv-(V'-) = 2sin2(V'_/2) = (1 
cosV'-), and HA.11() yields 



(COS(01 - 02)) = 1 



(|p«|2-l)(|p^|2-3) 



(A.13) 
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